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ABSTRACT 

' We present imaging and spectro- polarimetric observations of the ultraluminous 

^ ■ infrared galaxy IRAS P09104+4109 using the Keck 10-m Telescope. We detect the 

CN ' clear presence of broad H/3, H7, and Mg II A2800 emission lines in the polarized flux 

^ . spectra of the nucleus and of an extranuclear emission region ~ 4" away, confirming 

, the presence of a hidden central quasar. The polarization of the broad Mg II emission 

I line is high (~ 29%), consistent with the remarkably high polarization (~ 30%-40%) 

CS| I observed in the extended continuum emission. This indicates that the off-nuclear 

^ ' continuum is dominated by light scattered from the hidden quasar, most probably 

^ ■ by dust mixed with the line emitting gas. The high polarizations, combined with 

■ the "foreshortened" morphology of the polarized brightness distribution allow us to 

^f) . constrain the scattering biconical structure to be at inclination i ~ 50° with a 

. half-opening cone angle 9c ~ 40°. 

', The narrow emission lines are polarized in a stratified fashion, with the 

^ ; high ionization hues ([O III], [Ne V], [Fe VII]) being polarized 0.7%-1.7% and 

' [O II] essentially unpolarized. The line polarizations are positively correlated with 
1' 

I ' critical density, ionization potential, and velocity width of the emission lines. This 

2 ' indicates that, as is the case with the narrow-line radio galaxies, which also often 

c/3 . contain powerful quasars, the narrow-emission line region may be partially shadowed 

■ 

• • , by the putative torus, with the higher ionization lines originating closer to the nucleus. 

^ , 

. , One notable characteristic of the extranuclear knot is that all species of Fe are 

^ I markedly absent in its spectrum, while they appear prominently in the nucleus. In 

' addition, narrow Mg II is observed to be much weaker than predicted by ionization 

models. Our favored interpretation is that there is a large amount of dust in the 
extranuclear regions, allowing gaseous refractory metals to deposit. Near the nucleus, 
dust is destroyed in the strong radiation field of the quasar, inhibiting metal depletion 
onto grains. The extended emission regions are most likely material shredded from 
nearby cluster members and not gas condensed from the cooling flow or expelled from 
the obscured quasar. 
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The higher temperature inferred from [O III] hnes compared to that from [N II] and 
the general better agreement with models of line ratios, especially [O III] A5007/A4363 
and He II/H/?, provide strong evidence for matter-bounded clouds in addition to 
ionization-bounded clouds in the NLR. Ionization by pure velocity shocks can be ruled 
out. Shocks with photoionizing precursors may be present, but are probably not a 
dominant contributor to the energy input. 



Subject headings: galaxies: active — galaxies: individual (IRAS P09104+4109) — 
galaxies: Seyferts — infrared: galaxies — polarization 



1. Introduction 



Discovered by the IRAS all-sky survey, ultraluminous infrared galaxies (ULIRGs, see review 



by panders &: Mirabel 1996 ) emit most of their energy in the infrared {Ljji > 10^^ ^q)) s-'^d 



could harbor infant quasars enshrouded in a large amount of dust ( ganders et al. 1988[ ). On the 
other hand, they may also represent energetic, compact starbursts ( Pondon et al. 1991 ; Genzel et 



al. 1998| ). Much recent research effort has been devoted to understanding the dominant energy 
source in these ULIRGs - whether it is obscured quasars or intense bursts of star formation. 



IRAS P09104+4109 {z = 0.44; [Kleinmann et al. 1988]) is one of the few exceptionally 



luminous ULIRGs sometimes referred to as the "hyperluminous" [Ljji ^ 10 Lq) infrared 



galaxies. The others include F15307+3252 {z = 0.926; |Cutri et al. 1994| ) and F10214+4724 {z = 
2.286; Rowan- Robinson et al. 19911) . All of these galaxies have been shown to harbor a quasar 



nucleus obscured from direct view (P09104+4109, [Hines et al. 1999| F15307+3252, |Hines et 
al. 1995| ; F10214+4724, poodrich et al. 1996] ). This has led to the suggestion that perhaps all 
warm|^ ULIRGs contain buried QSOs and that they may be the misdirected type 2 QSOs (AGN 
with Seyfert 2-like emission-line characteristics but QSO-like luminosities). A diagnostic diagram 
(fTran et al. l999| ) involving [O III] emission line and infrared color /25//6O does indeed show that 
essentially all warm ULIRGs with sufficiently high ionization harbor energetic quasars in their 
centers. 



IRAS P09104-I-4109 has been identified with a central cD galaxy in a rich cluster ( [Kleinmann 



et al. 1988 ; Hall et al. 1997). It has a high-ionization spectrum characteristic of a Seyfert 2 galaxy. 
Ground-based and Hubble Space Telescope (HST) images show that IRAS P09104+4109 has an 
off- nuclear northern extension or "plume" of ionized gas ( [Kleinmann et al. 1988 ; Hutchings fc Neff 



1988 ; Armus et al. 1999 ) that may serve as a scattering mirror of the light originating from the 



^"Warm" ULIRGs are those having /25//6O > 0.2, |Low et al. 1988 ; 
IRAS flux densities in units of Jy at 25 /im and 60 /im, respectively. 



Sanders et al. 1988); f25 and /go are the 
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nucleus. We wish to probe the origin of this extension and constrain the scattering geometry for 
IRAS P09104+4109 by obtaining polarimetric observations of this knot. 

In this paper, we confirm the high polarization and the broad emission lines of II/3, II7 and 
Mg II in polarized flux of the nucleus. In addition to the S-band imaging polarimetry, we also 
present new spectropolarimetric data for the nucleus and off-nuclear emission regions of IRAS 
P09104-I-4109. We study in detail the emission-line spectra of the nuclear and extra-nuclear 
regions in order to derive an understanding of the ionization mechanism of the narrow-line region 
(NLR) of this galaxy. Throughout this paper, we assume Hq = 75 km s^^ Mpc~^, qo = and A 
= 0. At the redshift z = 0.44 of P09104-I-4109, 1" corresponds to a projected size of 5 kpc. 



Observations and Reductions 



Spectropolarimetric observations were made with the polarimeter ( pohen et al. 1997 ) installed 



in the low resolution imaging spectrometer (LRIS, pke et al. 1995 ) on the 10-m Keck I telescope 



on the night of 1994 December 31 (UT). A 1" wide, long slit, was centered on the nucleus of 
IRAS P09104-F4109, and oriented at PA = 14° to include the extended [O III] structure. The 
emission-line regions extend ~ 5" north- northeast of the nucleus ( Armus et al. 1999| ), as well 



as a few arcsecs south. We used a 300 grooves mm~^ grating which provided a dispersion of 
2.46 A pixel"^ and a resolution of ~ 10 A (FWHM). The observations were made by following 
standard procedures of rotating the half waveplate to four position angles (0°, 22.5°, 45°, and 
67.5°), and dividing the exposure times equally among them. Two sets of observations, lasting 60 
and 52 min., were made under photometric conditions, giving a total exposure time of 112 min. 
The results presented here are the average of the two sets. 

Spectra of the nucleus were extracted from an aperture 3 "4 wide, and those of the NE 
off-nuclear region were extracted from a region 3''6 wide centered 3^'4 from the nucleus. This 
includes the entire NE plume seen in the HST images by Armus et al. (1999). The data were 
reduced using VISTA, following standard polarimetric reduction techniques as described by Cohen 
et al. (1997). Hereafter, we shall refer to the nucleus spectrum as NUC and that of the NE 
extension as EXT. 

Imaging polarimetry was obtained on the night of 1995 December 16 with the Keck I 
telescope. The images were taken through the B filter, which contains mainly continuum emission. 
Exposure times were 10 min. at each of the four waveplate positions. Sky was clear but seeing 
was variable between 1"- 1''8. Reductions were done analogously to spectroscopic observations as 
described by Cohen et al. (1997). 



3. Results 
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3.1. Emission Line Ratios 

Table 1 presents the integrated emission-line flux ratios with respect to H/3 and their 
rest-wavelength equivalent widths for the NUC and EXT spectra, shown in Figure The emission 
line fluxes of NUC were measured from the starlight-subtracted spectrum. The contribution from 
an old stellar population represented by the elliptical galaxy NGC 821 was used to remove the 
starlight. The starlight fraction ranges from 5% in the blue end to a maximum of 40% in the red. 
We did not correct for starlight in the EXT spectrum as the galaxy contribution in the outer 
regions is negligible. The largest source of uncertainty in the flux measurement is the placement 
of the continuum. For strong lines (F/F(H/J) ^ 20%) the uncertainty is ^ 10%; for weaker lines 
(F/F(H/3) < 20%) the uncertainty is 10%-25%. 

The line ratios in NUC have been corrected for a small amount of extinction E(B — V)= 
0.24 determined from the H7/H/?, and H5/H/3 ratios, assuming the intrinsic case B values of 0.47 
and 0.26, respectively (Osterbrock 1989), and the extinction curve of Cardelli, Clayton, & Mathis 
(1989) with Ry = 3.1. This correction is small and does not significantly alter the conclusions of 
this paper. No corrections have been applied to EXT as its observed Balmer line ratios are close 
to case B values, indicating little extinction due to foreground dust. 

The small extinction derived might seem surprising, given the high infrared luminosity of this 
object. However, this is consistent with the result reported by Veilleux, Kim & Sanders (1999) in 
their spectroscopic study of a large sample of ULIRGs, that there is no obvious correlation between 
extinction and IR luminosity among the IR-bright galaxies. The mean E{B — V) for the Seyfert 2 
galaxies in their sample is 1.2 with a range of 0.3 to 2.7. Our reddening estimate is intermediate 
between those found by Kleinmann et al. (1988, E{B — V)< 0.12), who used the same optical line 
ratios, and Soifer et al. (1996, E{B — V) = 0.38-0.79), who used infrared [S II] and hydrogen lines 
in combination with optical lines. Soifer et al. also suggested that the low reddening may result 
from the dust being confined in a disk-like structure, which effectively obscures the central nucleus 
but does not block our line of sight to the NLR. As Veilleux et al. (1999) noted, the color excess 
derived from the optical line ratio method tends to underestimate the amount of dust in these 
objects. 

In the nucleus, the narrow lines all display a strong blueward asymmetry, as reported by 
Kleinmann et al. (1988), Crawford & Vanderriest (1996) and Hines et al. (1999), and cannot be 
adequately fit with a single Gaussian profile. Two Gaussian profiles including a main, central 
component with (intrinsic) FWHM ~ 600 km s^^ and another component blueshifted 720-980 km 
s~^ with FWHM ~ 700 km s~^ are able to fit the observed profiles well. Our measured range of 
velocity shifts of the blue-shifted component is somewhat smaller than the 1250 km s"^ reported 
by Crawford & Vanderriest (1996), and more consistent with 950 km s~^ reported by Kleinmann 
et al. (1988). 

Aside from the smaller line widths in EXT, the most remarkable difference between the NUC 
and EXT spectra is the complete absence of any Fe species in EXT, as Figure |I| shows. 
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3.2. Imaging Polarimetry 

The imaging polarization of IRAS P09104+4109 is shown in Figure |2|. The polarization is 
very high (~ 25%-30%) throughout the nucleus and extensions up to 5" away. The polarization 
increases strongly across the source. In the south, a number of pixels shows P > 30%, and the 
highest is 40.3% it 12.8%. We note that the distribution of surviving vectors is biased because 
those which have positive errors are preferentially saved; however, there is no doubt that the 
polarization is remarkably high. 

Unlike most of the narrow- line radio galaxies like 3C 195 and 3C 33 ( pohen et al. 1999| ), 
there is no well-defined V-shape bi-cone of scattered light, although the entire region is highly 
polarized. The polarization field displays a strong similarity to the patterns in the BLRG FSC 
2217+259 and the highly-polarized NLRG 3C 234 ( [Tran 1995| ; [Cohen et al. 1999| ). To the NE 



and SW, the vectors are perpendicular to the radius and this suggests scattering from a central 
illuminating source. However, the curvature of the polarization vectors is small, as in 3C 234 and 
FSC 2217+259. This may suggest that the viewing angle is small, indicating that the line of sight 
is just outside the ionization cone (see §4.3), a geometry that has also been suspected in 3C 234 
QTran et al. 1995| ; [Cohen et al. 19991 ). 



The southern extension appears to display somewhat higher polarization than its northern 
counterpart. This may be a result of forward scattering and would imply that the cone axis 



is pointed toward us in the south, analogous to Cygnus A ([Ogle et al. 1997[) . Dilution by line 
emission is not a major factor since there are no strong emission lines through the B band used. 
However, dilution by hot stars or nebular continuum could play an important role in the north as 
there is more ionized gas there. Our imaging polarimetry confirms and extends the polarization 
map obtained with HST by Hines et al. (1999), which shows a refiection "cone" in the central 1" . 

The right panel of Figure Q shows the polarized flux P x F which is commonly called the 
"polarized light" . As can be seen, the polarized light distribution extends more to the north than 
to the south, reaching up to ~ 4" (20 kpc) to the northeast along PA = 19°. This coincides well 



with the extended [O HI] emission seen in the ground-based ( Kleinmann et al. 1988 ) and HST 
images ( [Armus et al. 1999[) . The polarized extension is thus better aligned along this emission-line 
structure (PA ~ 11°) than the radio axis {PA = 333°), a condition also exhibited by many high 



and low-redshift radio galaxies ( [Tran et al. 1998[ ; Cohen et al. 199"9| ). This indicates that EXT 



serves as a reflection region of nuclear light, and that it lies exposed to the ionizing radiation 
escaping through the ionization cones from the hidden quasar. The polarized flux distribution 
peaks at the total flux peak, suggesting that scattered light dominates the observed radiation. 



3.3. Spectropolarimetry 
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3.3.1. Nucleus 

The spectropolarimetry for NUC is displayed in Figures ^. High polarizations are seen, with 
magnitudes consistent with those seen in the imaging polarimetry (Fig. ^). The Keck data of 
the nucleus confirm the general results of Hines et al. (1993, 1999). The observed continuum 
polarization uncorrected for starlight ranges from ~ 8% at 6000 A (rest) to about 18% at 3000 
A(rest). The polarization PA is constant with wavelength at 97°. In addition to broad H/3, H7, 
and probable Fe II multiplets near 3200 A reported by Hines et al., we also detect broad Mg II 
in polarized flux. We measured a FWHM of 7540 it 700 km s~^ for the broad H/3 in polarized 
flux, somewhat smaller than that reported by Hines et al. (1999, 12,000 it 2500km s^^). This 
is comparable to the FWHM of 6000 it 700 km s^^ measured for H7. The resemblance of the 
nuclear P x Fx to a quasar spectrum also extends to its shape. 

The galaxy-corrected total flux spectrum Fx {a = —0.94, f^, oc z^") is redder than 
P X Fx (a = —0.46), which is similar in slope to the more core-dominated composite QSO 
spectrum of Baker &: Hunstead (1995, a = —0.5), or that of Francis et al. (1991, a = —0.32). 
Our measured spectral indices agree well with those of Hines et al. (1999). This suggests that 
there is a second "featureless continuum" FC2 in the total flux spectrum that makes it redder. 
This component is likely due to a combination of hot stars and nebular continuum plus reddened 
transmitted light from the hidden AGN (e.g., Cohen et al. 1999| ). Quasars viewed at high 



inclination angle have been shown to exhibit redder continuum due to dust extinction ( [Bakei 



1997] ). The presence of hot stars, perhaps arising in a starburst, is suggested by the broad feature 
underlying He II A4686. This feature cannot be a broad scattered He II component since it appears 
very strongly in the total flux spectrum, suggesting that the light is viewed directly, and it is not 
detected in the P x Fx spectrum (see Fig. |^. It can be attributed to Wolf-Rayet (W-R) stars, as 
has been identified by Heckman et al. (1997) in Mrk 477 and Storchi-Bergmann, Cid Fernandes, 
k Schmitt (1998) in Mrk 1210, and Tran et al. (1999) in TF J1736+1122. 

In addition, there is nebular continuum giving rise to the Balmer jump near 3650A. If this 
is the case we would expect the Balmer jump to be unpolarized. Is it? To answer this, we have 
constructed an "unpolarized spectrum" = F^- (1/0.20) • P x Fx, assuming an intrinsic continuum 
polarization of 20%. We then measured the Balmer jump index defined as /3727- //3727+ , where 
/3727_ is the continuum flux just blueward of [O II] A3727, and 73727+ that just redward of it. We 
obtain an index of 1.74 it 0.1 for the total flux spectrum, and 2.0 ib 1 for the unpolarized flux. 
Naturally, there is a larger uncertainty in the value obtained from the unpolarized flux spectrum. 
However, a Balmer discontinuity appears to still be present, meaning that it is unpolarized, and 
supporting the notion that FC2 contains nebular continuum plus starlight from a starburst. 

We find that Ps, the intrinsic polarization^ produced by scattering, is Ps(H/3) ~ 16% while 
Ps{M.g II) ~ 29% (see below). This strengthens the reality of the nebular continuum discussed 



"We refer to the continuum-subtracted polarization of the broad emission Une as the "intrinsic polarization" 
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above, since a rising Pg toward the blue would leave most of the Balmer discontinuity in the 
unpolarized flux spectrum. This blue Pg is possible with dust but not electron scattering. Thus, 
dust scattering is important, while electron scattering could contribute but not dominantly. 

3.3.2. Extension Polarization Spectrum 

The polarization of EXT is noisier than NUC, and we have binned the results (except for 
F\), presented in Figure ^, by 5 pixels to improve S/N. As can be seen, the polarization shows 
similar behavior and magnitude to those in NUC. In addition, a broad Mg II is clearly seen in the 
total flux spectrum. This indicates that the continuum in EXT is dominated by scattered light 
from the obscured quasar. The polarization PA is roughly constant with wavelength at 106° ± 
2°, slightly higher than that measured in the nucleus. The PA is accurately perpendicular to the 
slit, as expected for an off-nuclear measurement of a reflection nebula (see [Cohen et al. 1999| for 
similar cases). 

We show the EXT P x Fx smoothed by 7 pixels in Figure |5|, which shows clearly that broad 
II/3 and Mg II are present. In P x F\, the spectral index is a = +0.5, significantly bluer than 
the average quasar spectrum. Assuming that the incident spectrum is typical for a quasar, this 
suggests that, as for the nucleus component, dust is the dominant scatterers which could bluen 
the light. 

3.3.3. Modeling of the Broad Line Polarization 

We model the spectropolarimetric observations using the procedure described by Tran et 
al. (1997), to separate the various components of the emission lines and continuum. In a region 
surrounding the [O III]+H/3 complex, we perform a fit of the continuum and emission-line 
components in the total flux, Q x Fx, and U x Fx spectra. The broad II/3 is assumed to have the 
same profile as in the polarized flux spectrum. The narrow lines of II/3 and [O III] were fitted 
assuming that they are similarly polarized, and consist of one main core component and another 
blue-shifted component, both with Gaussian profiles. Our Keck data do not have sufficient spectral 
resolution to warrant the addition of more than two components in the [O III] profile, and we do 
not attempt to fit the line polarizations with the three narrow components identified by Hines et 
al. (1999). The results indicate that the continuum and broad II/5 are similarly polarized at 16% 
± 3% and PA = 98°. 

Similar modeling at the Mg II line shows that the polarization of the broad Mg II line, at 
29%, is substantially higher than broad II/3, but is consistent with values seen in the imaging 
polarimetry. If real, it suggests that the intrinsic polarization rises to the blue, which is only 
possible if dust scattering is present. 
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Our fitting also shows that the core component of the narrow [O III] hnes has a small amomit 
of polarization, 0.70% it 0.03% at PA = 90°, which is different from that of the broad line and 
continuum. On the other hand, the [O II] A3727 line is virtually unpolarized. We measure a 
polarization of 0.22% it 0.12%. This indicates radial stratification of the NLR gas (see below). 



3.3.4- Stratification of Narrow- Line Polarization 

The high S/N data allowed us to measure the polarization of various emission lines in the 
nucleus. The results are in Table 2, where we list the (biased) P, 9, their associated errors, the 
ionization potential (IP) of the lower and upper stages of ionization for the ions of interest, and 
the critical density {ucrit) for collisional de-excitation of each line (calculated at 10^ K). The 
polarizations are flux-weighted means in the core component only, and the uncertainties on P and 
6 are la statistical errors alone and do not include systematics such as those due to continuum 
placement and data reduction. In Figure we plot the measured line polarizations versus Ucrit, 
the mean of the lower and upper IP, and FWHM. A clear trend is apparent in all three plots: 
lines of higher Ucrit, IP and FWHM exhibit higher polarization. Such a correlation has also 
been observed in the LINER NGC 4258 ( [Barth et al."l999| ). To our knowledge, NGC 4258 and 
IRAS P09104+4109 are the only two AGNs where this phenomenon has been shown to occur 
systematically over a number of forbidden lines with a wide range of ionizations. The tendency 
for [O III] to be more polarized than [O II], especially in radio galaxies ( |di Serego Alighieri et 



al. 1997 ), has been known for some time, and an anisotropic [O III] emission structure has been 



suspected as the cause ( Hes, Barthel, fc Fosbury 1993| ). Together with the well-known correlation 



of linewidth with Ucrit and IP (e.g., De Robertis &: Osterbrock 1986| ) in Seyfert galaxies, this 



result provides strong indication that the narrow emission-line structure in IRAS P09104+4109 is 
radially stratified, with lines of higher ionization and critical density being emitted closer to the 
nucleus, and therefore more likely to be partially occulted. The higher polarization results from 
less dilution by direct light. 



3.4. Physical Conditions and Diagnostic Diagrams 

The temperature determined from the [O III] A5007/[O III] A4363 ratio (using Table 1) is ~ 
1.3x10"^ K for both NUC and EXT. We also used the [N II] A6583 line from the near-IR spectrum 
of Evans et al. (1998) to measure the temperature from the [N II] A5755/A6583 ratio. This is 
only available for NUC, giving Te = 9000 K. This is about 4000 K cooler than the [O III] region, 
and may suggest that matter-bounded clouds are present (e.g., Wilson et al. 1997| ). The electron 



density estimated from the [Ar IV] A4712/A4740 ratio is ^ 10^ cm"^ and 5,000 cm"^ for NUC and 
EXT, respectively. This can be compared to a value of 6,000 cm~^ for NUC determined from [CI 
III] A5518/A5538 line ratio. The [Ar IV] AA4712/A4740 ratio is clearly better determined than the 
[CI III] A5518/A5538 ratio since the stronger Ar lines are less affected by starlight subtraction and 
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yield more certain flux measurements. 

To help us understand the rich emission-line spectrum of IRAS P09104+4109, we compare 
the line ratios with model calculations. In particular, we would like to understand the source of 
the ionization and the excitation mechanism of a large number of emission lines encompassing a 
wide range of ionization, from the low-excitation [N I] A5200, and [S II] to the coronal lines such 
as [Ne V] and [Fe VII]. Theoretical models from three main mechanisms are considered: 

1) Simple photoionization of essentially radiation-bounded NLR clouds by the central source 
having a powerlaw radiation spectrum. The emergent emission-line spectrum is characterized by 
the slope a of the ionizing continuum, and the ionization parameter U = Q/Anr'^cne, where Q is 
the number of ionizing photons, r is the radius and Ue is the electron density. We build the AGN 



models with the photoionization code MAPPING Ic ( Ferruit et al. 1997 ), using a power-law with 
index a = —1.5. The models are isobaric, so the density varies across the nebula. Recently, a new 
set of models has been computed for the NLR in an attempt to account for more realistic physical 
conditions. In these models the line emission originates from a wide range of temperatures and 
densities, although the observed emission-line spectrum results only from clouds best able to emit 
it. These are called "locally optimally emitting clouds" (LOG) models (Ferguson et al. 1997a, 
1997b). We shall discuss these models but compare our observations only with the simple, isobaric 
calculations that we made. 

2) Two-component photoionization model of Binette et al. (1996, 1997), consisting of both 
optically thin gas (matter-bounded, MB) clouds and optically thick gas (ionization-bounded, IB) 
clouds. The sequence is parameterized by the parameter ^m/i which is the ratio of the solid 
angle subtended by MB clouds relative to that of IB clouds, as viewed by the observer. Most of 
the higher ionization lines (i.e., coronal lines and [O III]) are emitted by MB clouds, which lie 
closer to the nucleus, filtering the radiation field that reaches the IB clouds further away, where 
essentially all of the lower ionization lines (e.g, [N II], [O I]) originate (Binette et al. 1996; Binette 



et al. 1997). A large value of ^m/i indicates a larger weight to the MB component, and therefore 



a higher excitation spectrum. 

3) Shocks and shock + precursors models of Dopita & Sutherland (1995, 1996). Mechanical 
energy deposited by velocity shocks due to turbulent cloud motions or the outflowing radio jets 
can produce powerful local UV radiation which can ionize the gas. Two types of shock models 
have been considered: pure shocks, in which mechanical energy from shocks is solely responsible 
for collisionally heating the gas, and the hybrid shocks + precursors, also known as "photoionizing 
shocks", in which the shock itself generates an abundant supply of ionizing photons in the hot 
post-shock gas and produce an extensive precursor H II region. 
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3.4- 1- Comparison to Model Predictions 

In Figures ^-|^, we present several diagnostic diagrams involving a number of emission lines 
with a wide range of ionization as an attempt to identify the main underlying mechanism for 
ionizing the NLR gas in NUC and EXT. Model calculations based on the above three mechanisms 
are shown for comparison with the observations. Simple power-law U sequences are plotted for 
density = 10^ cm~^ (of the front layer) as determined from the nuclear gas. Sequences with 
Uf. = 5000 cm~^ (as determined from the extended gas) produce very similar results and we do not 
plot them on the diagrams. Two different sequences are shown for the IB+MB photoionization 
models of Binette et al. (1997), each characterized by the value of the ionization parameter at 
the irradiated face of the MB component, Uo- high Uq = 0.5 and low Uq = 0.05. The ionization 
parameter of the IB component is 6.5 x 10^^ for both. On all figures, the parameter Ajy.j^j goes 
from 0.01 to 100, increasing from left to right. Pure shocks and shocks + precursors models are 
plotted as grids parameterized by the shock velocity V and the magnetic parameter B/y/n. V 
varies from 150 to 500 km s~^, and B/y/n varies from 1 to 4 /iG cm~^/^. All models are shown 
for solar abundances. 

The symbols have the following meanings: star denotes EXT and open circle represents 
the total line flux from both the core and blue-shifted components in NUC. The solid and open 
triangles indicate the core and blue-shifted emission line components of NUC, respectively. We 
shall now discuss these diagrams. 



3.4-2. Shocks and A GN photoionization 

Figure |^ shows that in two of the best optical diagrams for shock diagnostics, 
([O III] A5007/A4363 vs. He II A4686/H/3 proposed by [Clark et al.T997| , |Villar-MartiTrit 



|al. 19991 , and [Ne V]A3426/[Ne III]A3869 vs. [O III] A5007/H/? proposed by [Allen et al. 1998|) , 



the positions of the NUC component lie well outside the range of model grids predicted by pure 
shocks and shocks + precursors. The shock + precursor models are a little more consistent with 
the points for the EXT gas. These diagrams clearly rejects pure shock as the dominant emission 
line mechanism in both the nuclear and extra-nuclear regions of IRAS P09104+4109, as it tends 
to predict too low [O III] A5007. On the other hand, ionization-bounded AGN photoionization 
models have the opposite problem, failing to agree with the observed data by predicting too high 
[O III] A5007. 

The presence of dust (see §3.3, 4.1) in the EXT also weakens considerably any arguments 
that shocks play a dominant role in the main energy input mechanism. If shocks were present the 
dust grains would likely be destroyed (Donahue & Voit 1993; de Young 1999). Another argument 
against shocks is that the radio and emission line "cone" in IRAS P09104+4109 are misaligned 
by as much as 38° (Hines et al. 1999), making shocks induced by the nuclear jets unlikely. Also, 
the [O III] A5007/[O II] A3727 line ratio (a good indicator of the ionization parameter U) is 
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observed to be 5.4 for NUC and 2.11 for EXT. Thus U does not go up or remain constant in the 
extra-nuclear regions, as would be the case if shocks + precursors provide extra local heating (e.g., 
as in 3C 299, Feinstein et al. 1999| ), but actually drops, as expected from geometric dilution of the 



radiation field. Finally, the velocity widths of the emission lines in the extension are essentially 
unresolved < 550 km s~^, indicating that perturbations in the gas velocity by shocks, if any, are 
small. 

In conclusion, we can rule out pure shocks as the ionizing mechanism. Shock + precursors are 
also ruled out in the nuclear region, but it is more ambiguous whether they exist in the extension. 
Crawford & Vanderriest (1996) have suggested that shock heating may increasingly become 
important in the extended gas. However, we do not find sufficient evidence in our data to support 
a dominant role by shocks even in the extensions. More conclusive resolution between shocks and 
photoionization may come from additional diagnostics using UV lines ( [Dopita et al. 1997] ; [Allen et| 



al. 1998), which are more sensitive to shocks than optical lines. 



3.4-3. The Need for Matter- Bounded Clouds 

As Figure |^ indicates, the models for the high-density IB photoionization sequence predict 
too high [O III] A5007, typically giving [O IIIJ/II/3 ^ 16. This has been a well-known problem 
for IB photoionization models, known as the "temperature problem": the models predict too low 



a [O IIIJA4363/A5007 ratio and thus too low electronic temperatures (i.e., Binette et al. 1996| ). 
One attractive possibility that may solve this discrepancy is with the addition of MB clouds. 
As already mentioned, the presence of MB clouds is suggested by the higher temperature in the 
[O III] emitting region compared to that in the [N II] emitting region in the nucleus, since the 
MB clouds, where most of the [O III] originates, are expected to lie closer to the ionizing source 
QBinette et al. 1996| ; [Wilson et al. 1997] ). 



In Figure we present additional diagnostic line ratios with comparison to the IB -|- MB as 
well as shock -|- precursor models. The top two panels are the same as Figure 0. As can be seen, 
the shock -|- precursor models clearly fail to fit all the NUC data points and miss the EXT points in 
most diagrams. In general, the addition of MB clouds is able to resolve the two main difficulties of 
pure IB photoionization models: lines ratios of He II/H/3 and [O III] A5007/A4363. However, these 
sequences also have problems. For example, a common problem for all photoionization models is 
that [Ne V] is overpredicted (Fig. ^, top right). Note also that only models with A]\j/j < 1 are able 
to fit the observed He II, but they have difficulty reproducing the high [O III] A5007/H/3, especially 
the low Uo sequence (solid curve). The high Uq sequence (dotted curve in the plots) is consistent 
with both the observed He II and [O III] ratios, but it cannot reproduce the low-ionization lines [O 

II] /H/3, [N I]/H/3 and [S II]/H/3, which are predicted to be too strong. The low Uq sequence has 
the opposite problem, agreeing well with these low-ionization lines but failing to explain higher 
ionization lines like [Ne V] and [Fe VII], as well as the weakness of He II for the observed [O 

III] /H/3. Since it is the contribution from MB clouds that makes it possible to solve the He II and 
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temperature problems, and since the low-ionization lines are mostly due to the IB component, a 
higher ionization parameter for the IB component might be consistent with all data (recall that 



Binette et al. 1996 used a fixed U[IB]). 



In summary. Figure ^ shows that shocks + precursors models fails to reproduce the observed 
data points in NUC and, to a lesser degree, in EXT. Photoionization models with contribution from 
matter-bounded clouds are able to reproduced the data well, although no single photoionization 
sequence is perfect. 



4. Discussion 

4.1. Dust in the Extended Emission-Line Gas 

In §3.3, we presented polarimetric evidence for dust in the scattering regions. Diagnostic line 
ratios also show evidence for dust in the emission line gas. 



4-. 1.1. High- Ionization and Fe lines 

One of the most notable features of the emission line spectra of IRAS P09104+4109 is the 
conspicuous absence of all species of Fe in EXT and their strong presence in NUC (see Fig. |l| and 
Table 1). No other atomic species other than iron is seen to exhibit this behavior. This dramatic 
difference in Fe line strength between NUC and EXT cannot be due to the lack of high ionization 
photons in EXT, for we see clearly in EXT strong high- ionization lines of [Ne V], which has 
comparable ionization potential (97.1 eV) to [Fe VII] (100.0 eV). Furthermore, lower ionization 
Fe lines like [Fe VI], [Fe V] and perhaps [Fe III] are all completely absent from EXT. Coronal 
lines can be expected to be seen in the tenuous off- nuclear extended gas ([Ferguson et al. 1997b| ), 
and indeed [Ne V] is seen in our EXT spectrum, and [Fe IX] has been observed in the extended 
emission line regions of Circinus (|01iva et al. 1994 ; [Moorwood et al. 1996| ), Cygnus A ( pgle et 



al. 19971) and Tololo 0109-383 (Murayama, Taniguchi, k Iwasawa 1998). 



The peak range of the line strength distribution function in the LOC models of Ferguson et 
al. for [Ne V] and [Fe VII] overlap substantially with each other. Therefore, under "ordinary" 
NLR conditions (i.e., solar abundances, free of dust) we would expect to see [Fe VII] wherever 
[Ne V] is seen. If [Fe VII] is as strong relative to other lines as observed in the nucleus, there 
should be no problem detecting it in EXT. In NUC, [Fe VII] A6087 is about 3 times weaker than 
[Ne V] A3246 and 5.6 weaker than II/5. At this level, [Fe VII] should be even stronger than [N I] 
A5200, or about as strong as He I A4471, both of which are clearly seen in EXT. Even if it is 6-13 
times weaker, as predicted by models of Binette et al. (1997) or Ferguson et al. (1997), it should 
be at least about as strong as [Ar IV] AA4712, 4740 or He I A5876. Somehow, virtually all gaseous 
Fe is removed from the extension gas. With the inclusion of dust grains, [Fe VII] is expected to be 
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about 74 times weaker than [Ne V] (Ferguson et al. 1997b). At this level, it would be impossible 



to detect with our data. This strongly suggests that there is significant amount of dust in the 
extension allowing gas-phase Fe to deplete onto grains. This notion has been invoked to explain 
the emission-line spectra of Seyfert galaxies (Kingdon, Ferland, & Feibelman 1995). Neon is an 
inert gas and not expected to be depleted (see model calculations by Ferguson et al. 1997b| ). 



[Fe VII] A6087 is a coronal line, and both Ferguson et al. (1997b) and Binette et al. (1997) have 
suggested that coronal lines are not efficiently emitted in gas with dust. The lack of [Fe VII] A6078 
in EXT therefore is consistent with the presence of dust. 

It is unlikely that the absence of Fe lines in EXT is due to the lower metallicity of the EXT 
region, for other metals such as oxygen and nitrogen seem to be well represented. We would 
expect Fe to show a normal abundance as there exists a well-known [Fe/H] to [O/H] relation for 
the nebula gas observed in the solar neighborhood and galactic bulge (e.g., Maciel 1999). 



4.1.2. Narrow Mg II 

Another element which may lend itself well to diagnostics for dust is Mg since it is about as 
sensitive as Fe to depletion. Unlike the higher-ionization forbidden Fe lines, Mg II is complicated 
by several factors. It is a resonance line so it can be strongly affected by intervening dust 
extinction and gas absorption. On the other hand, like other low-ionization lines such as [S II], 
[O II], [N II] it may even get strengthened by photoelectronic heating of the gas by dust grains 
( [Ferguson et al. 1997i^ ). 

With these caveats, we now consider whether the resonant character of Mg II A2800 and its 
strong sensitivity to dust and gas absorption suggest the presence of dust as well. Figure ^ shows 
the diagram of Mg II A2800/[O II] A3727 vs. [O III] A5007/H/3. The MB + IB photoionization 
models predict Mg II to be stronger than observed in both the nuclear and extended gas. Shock 
-|- precursor models show reasonable agreement to the data, but they fail for most other lines, and 
have been eliminated for the nuclear emission (see §3.4.2, 3.4.3). One way to weaken Mg II A2800 
without dust is with intervening absorbing gas, containing Mg II. Since Mg II is in ionized form, 
the most likely explanation is that the absorbing gas is in the ionized cones. We compare the 
predictions of the photoionization models taking into account viewing effects: depending on the 
viewing angle of the observer relative to the clouds, the spectrum of the resonant lines, like Mg II 
A2800 will be different (Villar-Martin, Binette & Fosbury 1996). If the clouds are viewed from the 
illuminated face, the Mg II emission will be relatively stronger than when seen from behind. In 
the former case, the photons find a lower column density of Mg II (since the gas is highly ionized 
towards the illuminated face) and they can escape the nebula freely. However, photons that try 
to leave the cloud through the "non-illuminated" face of the cloud will encounter a much larger 
column density of Mg II and will be repeatedly absorbed. 

In Figure ^ we show this effect. The upper thick solid line shows the line ratios for a nebula 
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viewed from the illuminated surface. The lower thick solid line shows the line ratios for a nebula 
seen from behind. As can be seen, the difference is large. These models show that we need to view 
the clouds directly from the back (the non-illuminated side) to reduce the Mg II emission to the 
low observed values. However, this extreme view is probably not the case, for if the clouds are 
seen from behind, the observer would be inside the ionization cone and therefore, would have a 
full view of a quasar or broad line radio galaxy. Another possibility is that the absorbing gas is 
outside the ionized cones, but in this case the ionization of the gas cannot be explained. 

Thus again, the most reasonable explanation is that the emitting gas contains dust. It could 
be due to depletion (Mg is about as sensitive to depletion as Fe), or to a combination of resonant 
scattering and dust absorption. Pure resonant scattering alone cannot explain the weakness of the 
line. This conclusion is consistent with the indication for dust from the absence of Fe lines, as well 
as the polarized flux spectrum, which is blue compared to a composite quasar spectrum (§3.3.2). 

4.2. Origin of the Extended Gas 

While the indication for dust in the nucleus and extended emission regions of IRAS 
P09104+4109 comes as no surprise from one of the most luminous infrared emitters, this 
confirmation shows that the existence of dust is extensive over a large scale, and could have 
important implications for the origin of the extended gas. 

As suggested by Donahue k, Voit (1993) for the intracluster nebular gas that they studied, 
the presence of dust inferred in the extension implies that EXT cannot be gas condensed directly 
out of the X-ray emitting intra-cluster medium (ICM). This is because dust simply cannot survive 
in the hot ICM in any appreciable amount. Any dust introduced in the hot ICM is expected to 
be destroyed long before the gas is cooled enough to be visible or gets transported to the cluster 
center by the cooling flow. Villar-Martin &; Binette (1997) reached a similar conclusion in their 
flnding of dust in the extended emission-line regions of radio galaxies. 

Since there is no evidence that the metal abundances are anomalously low (§4.1.1), the most 
likely candidate for the extended gas is perhaps shreds of unfortunate members of the cluster 
being cannibalized by IRAS P09104-I-4109. Active star forming regions, either in the galaxies 
themselves or triggered by the interactions, are where copious amounts of dust are found. This 
in turn implies that the extranuclcar gas is unlikely to be material expelled or outflowing from 
the nucleus of the central quasar. HST WFPC2 imaging of IRAS P09104-F4109 by Armus et 
al. (1999) does indeed show what seem like whiskers and filaments, suggestive of stripped remains 
from nearby cluster members. In addition, recent HST WFPC2 and NICMOS imaging studies of 
ULIRGs by Surace et al. (1998) and Scoville et al. (2000) have shown that many galaxies contain 
numerous circumnuclear clusters of bright young stars (age ~ 10^-10^ yrs), whose forming regions 
are often associated with dust. 

The external origin for the gas is supported by the fact that many ULIRGs show evidence for 
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interaction or merging (e.g., Due, Mirabel & Maza 1997; Borne et al. 2000; gcoville et al. 200q ), 
which has been suggested to stimulate the ULIRG phase. Velocity structure of the extended 
nebula of IRAS P09104+4109 derived from integral field spectroscopy by Crawford & Vanderriest 
(1996) also suggests that the gas is probably external in origin. However, these authors do not 
favor the picture of the gas being merged with the cD galaxy, based on the low relative velocity 
(~ 100-200 km s~^) with respect to the nucleus. 



4.3. Scattering Geometry 

The polarizations of ULIRGs with hidden broad-line regions are generally very high 
(~ 20%-30%), suggesting that the inclination angle could be fairly substantial ( ^ 35°). Under the 
assumption that the lack of a well-defined V-shape polarization fan suggests that the inclination 
angle i does not exceed the half-opening cone angle 6c by a large amount, we can use the 
morphology of the polarization image (Fig. |2|) to constrain the scattering geometry of IRAS 
P09104-F4109, using the model of Brown & McLean (1977) (see also [Miller fc Goodrich 1990| ; 



Miller et al.lJ9T| ; [Wills et al. T99^ ; [Balsara fc Krolik 1993[ ; [Brotherton et al. 1998| ). Analysis by 



Mines et al. (1999) has constrained the viewing inclination and half opening cone angle lie in the 
range 34° < i < 41° and 15° < 6c < 33°, based on P= 20%. With our higher polarization of ~ 
30% from the broad Mg II and imaging polarimetry, these limits are still consistent, but we can 
eliminate the combination of both upper limits or lower limits for i and 6c (i.e., i = 41° 6c = 33°, 
or i = 34° dc = 15°), as they give too low P. The same applies for the combination i = 34°, 
6c = 33°, in which case our line of sight is just outside the cone, but P is only 15%. P agrees 
better with observations if i and 6c are closer to the upper and lower limits, respectively, but in 
this case their difference i — 6c = 29° is too large. If we limit the half opening cone angle to be 
^ 40° constrained by the HST image of Hines et al. and Armus et al. (1999), the difference i — 6c 
is minimized to ^ 10° for i ~ 50° and 6c ~ 40°. Therefore, this appears to be the best scattering 
geometry for IRAS P09104-I-4109 consistent with observations. 



5. Summary and Conclusions 

Our high-quality Keck spectropolarimetric and imaging polarimetric data confirm the high 
polarization and broad Balmer emission lines in the polarized flux spectra of IRAS P09104-I-4109. 
In addition, we also detect broad Mg II in both total and polarized flux spectra of the nuclear 
and extension regions, indicating the presence of a hidden quasar visible only in scattered light. 
The narrow-line polarizations exhibit a strong positive correlation with line width, critical density 
and ionization potential of the transition, indicating that the line emission arise from a radially 
stratified gas. The lack of a clear fan-like morphology in the polarization image suggests that our 
viewing angle (i 50°) is not far outside the ionization cone {6c ~ 40°) of IRAS P09104-I-4109. 
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The emission-line spectra of IRAS P09104+4109 are consistent with AGN photoionization 
with contribution from both ionization-bounded and matter-bounded clouds. Under the 
assumption that photoionization by the central hidden quasar is the main operating mechanism, 
[Fe VII] should be detected in the extended gas and Mg II should be much stronger than observed. 
The absence of all Fe lines and weakness of Mg II in the extension can be readily explained if 
dust is mixed with the ionized gas of the NLR. Furthermore, the fact that both the nuclear and 
the extended continua are highly polarized with relatively blue polarized flux spectra compared to 
that of the average quasar also supports the existence of dust, which could efficiently scatter and 
bluen the light. The other possibility is that the extended gas simply is highly deficient in iron, or 
is metal poor. This does not seem very likely, since the strengths of the other metal lines appear to 
be in good agreement with solar abundance. Since grains remain intact, this in turns implies that 
shocks are probably not important in the extra-nuclear regions. Although the diagnostic diagrams 
offer somewhat inconclusive assessment of the shocks+precursor models for the extension region, 
there is no strong evidence for shocks in IRAS P09104-I-4109 from the line ratios. Ionization by 
pure shocks can be ruled out by the data. The existence of dust and lack of strong shocks in the 
extended gas also suggest that it is probably dismembered remnants of a cluster neighbor being 
disrupted or in the process of merging with the central cD galaxy. 

Two predictions follow directly from our interpretation of the current data: The near-IR line 
[Fe II] 1.257 i^m has been observed in the nucleus of IRAS P09104+4109 by Soifer et al. (1996), 
but should be absent from the extension, confirming the absence of gaseous Fe, and hence the 
lack of shocks and presence of dust there. Likewise, [Ca II] AA7291, 7324 should be absent in 
EXT spectrum, since Ca is more than 10 times more sensitive to grain depletion than Fe or Mg 
([Ferguson et al. 1997b| ). 
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Fig. 1. — Comparison of the total flux spectra of the nucleus (top) and off-nuclear NE extension 
region {bottom). "X" denotes where Fe line is expected but not seen in the extension spectrum. 
Note the complete disappearance of Fe lines, in order from left to right: [Fe VII] A3588, [Fe V] 
A4229, [Fe VI+VII] AA4146, 5159, [Fe III+VII] AA5270, 5276, [Fe VI] AA5424, 5427, [Fe VII] A5720, 
[Fe VII] A6087. 
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FSC 09104+4109 




Fig. 2. — Polarization map for IRAS P09104+4109 obtained in S-band with the Keck telescope. 
The left panel shows the polarization vectors superimposed on the contour plot of the total flux, 

with north up and east to the left. The scale is 2" per tick, and the contours are 0.5, 1, 2, 4, 10, 
20, and 50% of the peak. Polarization vectors are shown for those pixels where P > 2.5(7. The 
line marked "R" shows the radio axis. A horizontal bar representing 100% polarization is shown 
at the bottom right. The right panel shows the polarized flux P x F. The contours are 5, 15, 
30, 50, and 80% of the peak. The polarized light distribution peaks at the flux peak, suggesting 
that the observed light is dominated by scattered radiation. The 1" wide spectroscopy slit ("S") is 
indicated, along with the extraction apertures for the NUC and EXT spectra. 
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Fig. 3. — Spectropolarimetry of the nucleus of IRAS P09104+4109. (a) Total flux spectrum Fx, 
displayed at two scales to show the weaker emission lines, (b) observed degree of polarization P, 
and (c) polarized flux spectrum P x F^. The flux scales are in units of 10~^^ ergs cm"^ A~^. 
Broad H/3, H7, and Mg II, as well as narrow lines of [O III], [Ne V], and [Ne III], are clearly present 
in the P x Fx spectrum. 
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Fig. 4. — Spcctropolarimetry of the NE off-nuclear emission line region of IRAS P09104-I-4109. 
From top to bottom panels are: the total flux, P, and polarization PA. P and PA have been binned 
5 pixels and the flux has been binned 2 pixels. The polarization shows similar high magnitude and 
wavelength dependence as in the nucleus. 
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Fig. 5. — The polarized flux spectrum P x F\ of the extension, smoothed 7 pixels to improve S/N. 
A broad H/3, Mg II and perhaps H7 are visible atop a blue continuum with a = +0.5 (solid line, 



-26- 



Oh 



2 

1.5 

1 

0.5 




Tmi| I I I Mlllj I I I IIMIj I I I MINI ^ I M I 



[N I] 



[0 111] 



viH — 



[0 II] 

I I I I 



[Ne V] 



[Fe 
[Ne III] 



I ' ' ' I 



I I I 



I I I I I I I I I I I I I I I I I 



I I I 



lO-* 105 IQ6 IQ7 20 40 60 80 100 120 
n j, (cm-3) mean IP (eV) 



2 - 



1.5 



0.5 



I ' ' ' ' I 



T 



I !^ 



I — 14 ■ I 

I — i — I 



I 1 1 



I I I I I I I I I I I I I I I I I I I I I I I I 





300 400 500 600 700 800 
FWHM (km s'l) 



Fig. 6. — Narrow forbidden line polarization as a function of ricrit, ionization potential (IP) and 
line widths. A clear trend of higher polarization with higher ricritt 

IP, and FWHM is evident. 
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Fig. 7. — Predictions of shock (dotted lines), shock+precursor (dashed Hnes) and AGN 
photoionization models with radiation bounded clouds and density Ue = 10^ cm~^ (solid lines) . U is 
the ionization parameter that defines the AGN sequence (see text). U varies between 0.001 and 0.46. 

Each shock (and shock+precursor) sequence is defined by a fixed value of the magnetic parameter 
B (B = B I ^fn) and a range of shock velocity v (150-500 km s^^). Three shock sequences are 
shown, corresponding to B = 1-4 //G cm~^/^. These diagrams clearly reject pure shocks (radiative 
cooling of shocked gas) as the dominant emission line mechanism in the nuclear and extended gas. 
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Fig. 8. — Shock+precursor model and AGN photoionization model (with contribution from 
matter bounded clouds) predictions. The top two diagrams are the same as in Figure |7[ The 
shock+precursor models are plotted only on those diagrams where the calculated line ratios are 
available. Symbols are as in Figure 0. The filled circles in the MB sequences represent the models 
with A]^f/j = 1 and the arrow in the upper panels shows increasing Aj^/j. The models with matter 
and ionization bounded clouds are in good agreement with most line ratios except [Ne V]/[Ne 
III], whose predicted values are too high compared with the measurements. These diagrams 
reject shock+precursor models for the nuclear emission, where, therefore, AGN photoionization 
dominates. The discrimination between AGN photoionization and shock+precursor models is not 
possible for the extended gas. 
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log([OIII]/Hb) 

Fig. 9.— Mg II A2800/[O II] A3727 vs. [O III] A5007/H/3. The thick sohd hnes show the 
photoionization U sequence with radiation bounded clouds as shown in Figure but taking into 
account perspective effects. The upper sequence corresponds to the clouds seen from the illuminated 
face and the lower sequence corresponds to the clouds seen from the rear. In this scenario, only 
the extreme case such that we are inside the ionization cone and see the clouds from the rear could 
explain the weakness of Mg II A2800. See text for details. For comparison, shock + precursor and 
MB photoionization models are also shown. Symbols and lines are as in Figure |^ 
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Table 1. Emission Line Flux Ratios, Equivalent Widths and FWHMs 



IRAS P09104+4109 Nucleus" IRAS P09104+4109 Extension 



Line Flux Ratio^ EW'^ FWHM'^ Flux Ratio'' EW= FWHM' 



He II A2733 

[Mg V] A2783 0, 

Mg Iln A2798 0, 

Mg lib A2798 0, 
He I A2830 + O III A2836 0, 

[Ar IV] A2854 0, 

[Ar IV] A2868 0. 

[Mg V] A2928 0, 

He I A2946 0, 

O m A3047 0, 

O HI A3133 0, 

He I A3188 0, 

He II A3203 0, 

O HI A3312 0, 

[Ne V] A3346 0, 

[Ne V] A3426 0, 

O HI A3444 0, 

[Fe VII] A3588 0, 

[O II] A3727 1, 

H12 (+ Hll) 0, 
Hll 

HIO 0, 

H9 

[Ne HI] A3869 0, 

H8 + He I A3889 0, 
He A3970 + [Ne III] A3967 

He I A4026 

[S II] A4071 0, 

R6 A4102 0, 

He I A4144 0, 

[Fe V] A4229 0, 

H7 A4340 0. 

[O III] A4363 0. 

He I A4471 0, 

He II A4686n 0, 

He H A4686b 0, 

[Ar IV] A4712 0, 

[Ar IV] A4740 0, 

H/3 A4861 1, 

[O III] A4959 3, 

[O III] A5007 1 
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5.93 




,0667 


9.25 




,8060 


121.7 


725, 465 
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27.9 
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,3824 


62.3 


628, 580 


,0112 


1.8 




,09695 


15.9 


712, 750 


,2163 


35.4 


835, 465 


,00980 


1.6 




,07158 


12.0 


918, 859 


,4127 


70.8 


731, 480 


,1588 


27.1 


567 


,0359 


6.4 


906, 433 


,1518 


27.6 


400 


,1861 


34.0 


5020 


,0150 


2.7 




,0193 


3.5 




,000 


184. 


660, 465 


,794 


1050. 


687, 495 


1.18 


3541. 


700, 465 



0.3110 


17.4 


840 


0.71053 


55.6 


14800 


0.0324 


1.84 




0.0416 


2.74 




0.0367 


2.79 




0.0361 


2.78 




0.06463 


5.09 




0.07033 


5.65 




0.1561 


12.7 




3.646 


377. 




0.0326 


3.82 




0.0408 


5.1 




0.0588 


8.08 




0.0734 


10.2 




0.8343 


123.3 




0.1918 


29.6 




0.4041 


59.2 




0.0203 


2.9 




0.1145 


17.0 




0.2705 


41.1 




0.0134 


2.05 




0.4655 


73.9 




0.1343 


22.0 




0.0417 


6.40 




0.1169 


17.4 




0.0245 


3.73 




0.025 


3.90 




1.000 


140.0 




2.574 


369.5 




7.701 


1084. 
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Table 1 — Continued 



Line 


IRAS P09104+4109 Nucleus" 


IRAS P09104+4109 Extension 


t lux Katio 


jhjVV r WrlM 


Flux Ratio'' EW^ FWEM"^ 


[Fe VI] A5146 


0.0318 


6 36 567 




[Fc VII] A5159 


0.0376 


8.0 560 




[Fe VI] A5176 


0.0289 


5.82 




[N I] A5200 


0.05011 


11.0 


0.1436 24.9 


[pe TTT -1- Fp VTT + Fe TTl^ 


0.0378 


8.3 




[Ca V] A5309 


0.0361 


8.0 




[Fe II + Fe VI] •'^ 


0.0225 


5.0 




[Fe III + Fe VI] » 


0.0436 


9.96 




[Fc VI] A5485 


0.0142 


3.26 




[CI III] A5518 


0.0117 


2.46 




[CI III] A5538 


0.0149 


3.14 




[Ca VII] + [Fe VI] 


0.0209 


4.72 




[Fe VI] A5677 


0.0153 


3.44 




[Fe VII] A5720 


0.1142 


26.5 




[N II] A5755 


0.0137 


3.73 


0.0359 5.91 


He I A5876 


0.1174 


29.5 680, 510 


0.1028 17.7 


[Fe VII] A6087 


0.1805 


46.3 774, 674 





'^Flux ratios and EWs are measured from starlight-subtracted spectrum. 

''Observed line flux ratios relative to narrow H/3. Values for the nucleus are the sum of both the central 
and blueshifted components corrected for E{B — V) = 0.24. 

■^Rest-frame equivalent widths in A. 

■^FWHMs (km s~^) have been corrected for instrumental resolution of 550 km s~^. The two values for 
the nucleus denote the central and blueshifted component, respectively. 

''Probable identification: [Fe III] A5270 + [Fe VII] A5276 + [Fe II] A5262. 
fProbable identification: [Fe II] A5334 -|- [Fe VI] A5335. 
eProbable identification: [Fe III] A5412 + [Fe VI] AA5424, 5427. 
''Probable identification: [Ca VII] A5615 + [Fe VI] A5631. 
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Tablc 2. Forbidden Line Polarization 



Line 


P(%) 


PA (°) 


Lower IP (eV) 


Upper IP (eV) 


Tlcrit (cm ^) 


[Ne V] A3346 


1.77 ± 0.38 


104.0 ± 6.2 


97.1 


126.2 


1.38 X lO'^ 


[Ne V] A3426 


1.06 ± 0.17 


102.4 ± 4.6 


97.1 


126.2 


1.38 X 10^ 


[0 II] A3727 


0.23 ± 0.12 


124.8 ± 15.2 


13.6 


35.1 


4.20 X 10^ 


[Ne III] A3869 


0.70 ± 0.07 


93.8 ± 3.0 


41.0 


63.5 


8.48 X 10^ 


[Fe V] A4229 


0.63 ± 0.43 


56.9 ± 17.7 


54.8 


75.5 




[0 III] A4959 


0.70 ± 0.03 


87.7 ± 1.1 


35.1 


54.9 


6.16 X 10^ 


[0 III] A5007 


0.63 ± 0.03 


85.8 ± 1.3 


35.1 


54.9 


6.16 X 10^ 


[N I] A5200 


0.90 ± 0.77 


118.0 ± 24.6 


0.0 


14.5 


2.34 X 10^ 


[Fe VII] A6087 


1.42 ± 0.71 


68.4 ± 14.4 


100.0 


128.3 


3.6 X 10'^ 



